Crossovers from parity conserving to directed percolation universality 
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The crossover behavior of various models exhibiting phase transition to absorbing phase with par- 
ity conserving class has been investigated by numerical simulations and cluster mean-field method. 
In case of models exhibiting Z2 symmetric absorbing phases (the NEKIMCA and Grassberger's 
A stochastic cellular automaton) the introduction of an external symmetry breaking field causes 
a crossover to kink parity conserving models characterized by dynamical scaling of the directed 
percolation (DP) and the crossover exponent: l/cp ~ 0.53(2). In case an even offspringed branching 
and annihilating random walk model (dual to NEKIMCA) the introduction of spontaneous particle 
decay destroys the parity conservation and results in a crossover to the DP class characterized by 
the crossover exponent: l/cp ~ 0.205(5). The two different kinds of crossover operators can't be 
mapped onto each other and the resulting models show a diversity within the DP universality class 
in one dimension. These 'sub-classes' differ in cluster scaling exponents. 
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I. INTRODUCTION 

The study of nonequilibrium phase transitions is an im- 
portant task of statistical physics. Genuinely nonequilib- 
rium transitions can be observed most easily in models 
exhibiting transition from an active to an "absorbing" 
state, where the fluctuations are negligible, hence no re- 
turn is possible. The exploration of critical phenomena 
and universality classes of simple, one component models 
has been started [E S H, 0] and important steps towards 
a full classification have been done 

For a long time it was a common belief that all con- 
tinuous, nonequilibrium phase transitions belong to class 
of the directed percolation (DP) and a hypothesis ad- 
vanced by Janssen and Grassberger 0, H, Q ■ This states 
that in one component systems exhibiting continuous 
phase transitions to a single absorbing state (without 
extra symmetry and inhomogeneity or disorder) short 
ranged interactions can generate DP class transition only. 
Despite the robustness of this class experimental obser- 
vation is rare, owing to the high sensitivity to disorder 
and long-range interactions. A very recent experimental 
study [lOj has reported clear and comprehensive experi- 
mental evidence of DP criticality. Later it was discovered 
that in systems with infinitely many frozen absorbing 
states (IM type) [HI [H, HI, [3 like in the pair contact 
process (PGP) the static exponents coincide with those 
of the DP. The dynamical cluster spreading behavior is 
different owing to the long time memory generated by 
the frozen monomers [l5| . In the literature the transi- 
tion type of PGP is often called DP type. Very recently 
(l6l | have investigated in 1 + 1 dimension the crossover 
from PGP to DP type of models by introducing differ- 
ent absorbing state reduction mechanisms and nontrivial 
exponents have been found. This finding confirms that 
strictly speaking the universality class of PGP and DP is 
different. 

The first example for clearly non-DP critical behav- 
ior was found among stochastic cellular automata (SGA) 



by Grassberger [T3| . These models exhibit Z2 symmetric 
absorbing states and an effective kink dynamics, which 
follows an even offspringed branching and annihilating 
random walk (BARWe). In reaction-diffusion (RD) par- 
ticle models, with BARWe dynamics {A 3A, 2A 0) 
such phase transition class was discovered by [H, Ell. 
Since this is different from odd offspringed branching and 
annihilating random walks, in which DP class transition 
occurs, the name "parity conserving" (PG) was intro- 
duced to denote this class ^. 

An important example of the PG class behavior was 
discovered in one-dimensional kinetic Ising models with 
combined zero temperature spin-flip and finite temper- 
ature spin-exchange dynamics (NEKIM) [20!|. Here the 
domain walls between up and down spins follow BARWe 
dynamics and an exact duality transformation in one di- 
mension between the NEKIM and the BARWe particle 
model was established by [lH 2. Naturally the NEKIM 
exhibits two, Z2 symmetric absorbing states. This uni- 
versality class has not been observed in nature yet, how- 
ever we pointed out [2^ that this model is very insensitive 
to the quenched disorder, therefore it is a good candi- 
date for experimental veriflcation. In fact it's inactive 
phase, where annihilating random walk dominates has 
been observed 23] unless very strong disorder fractures 
the medium. 

The introduction of a symmetry breaking external 
magnetic field, in NEKIM, which favors one of the ab- 
sorbing states (but preserves the BARWe dynamics) was 
shown to change the type of transition from PG to DP 
type [2^ . Such crossover mechanism has been observed 
in other Z2 symmetric models as well (H HI, [H [H, Hi] . 



^ However other different names for this class like "Directed 
Ising" (DI) or "Generalized voter (GV)" or "BARW" can also 
be found in the literature 

^ Note that such mapping is not possible in higher dimensions, but 
since the upper critical dimension is dc = 4/3 it is not important. 
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On the other hand simulations [30| and field theory 
[31I [3^ proved that a PC breaking ^ ^ reaction in 
the BARWe model should also change the type of phase 
transition from PC to DP type. 

The more detailed study of crossover behavior among 
nonequilibrium universality classes has been intensified 
in the recent years. The universal crossover exponent, 
defined by the shape of phase boundary {rc{w)) as one 
introduces a relevant scaling field (w) (see for example 



Tciw) 



,1/0 



(1) 



has been determined in case of DP to the compact DP 
class [IIjIsIjIH,!!^]- Crossover between DP and isotropic 
percolation was investigated by field theory [sl, 
and simulations [4l|. Very recently the numerical ex- 
ploration of the crossover behavior from the diffusive 
pair contact process (PCPD) to DP [i^l has strength- 
ened the existence of the independent PCPD class the- 
ory. In IM type models 'nontrivial crossovers to DP class' 
have been found [l^ • The crossover behavior from DP to 
mean-field, generated by long-range diffusion, following 
the early studies [H, |3| ■, has currently been re-examined 
via numerical techniques (45| and field theory [46] and es- 
timates for the exponent (j) have been provided. Similar, 
diffusion driven crossover in case of PC class was deter- 
mined within the framework of NEKIM [47] long ago. 

In the present paper we determine the crossover expo- 
nent from PC to DP class in a SCA version of the NEKIM 
model [11] (NEKIMCA) and compare it with other re- 
alizations of the PC class. In particular we confirm the 
universality of 4> in case of Z2 symmetry breaking fields 
by simulating Grassberger's A model. We compare this 
crossover behavior with the outcome of parity conserva- 
tion breaking in a BARWe model. 



II. CROSSOVER OF NEKIMCA IN AN 
EXTERNAL FIELD 

The NEKIM exhibiting PC class transition was sug- 
gested by ^2(| as a generalization of the Glauber Ising 
model [4g|. It is defined by the alternating application 
of a T = spin-flip sweep and a T > Kawasaki spin- 
exchange update of a one dimensional lattice. While the 
spin-flip dynamics generates annihilating random walk 
of kinks, the spin-exchange introduces a parity conserv- 
ing branching {A 3) of the domain walls. Tuning the 
relative strengths of the reactions one can get a phase 
transition from an active to a kink-free, adsorbing state 
(the order parameter is the density of kinks) . 

It was realized in f23| that branching reactions appear 
automatically in the SCA version of the NEKIM, if the 
spin-flip update is done synchronously due to overlaps. 
With this dynamics we can obtain a very simple model, 
with PC type of critically, which can be implemented on 
a computer by efficient bit-coding. The NEKIMCA spin 
updates, represented by bit field operations of a computer 



word, generating BARWe reactions of the kinks are the 
foUowings: 

• Random walk of domain walls (•): 

T T-i^T-i i 

generated by a spin-flip between oppositely ori- 
ented spins, with probability Wi, 

• Annihilation of a pair of kinks: 

T-i-T ^ T T T ^ 

generated by a spin-flip between a pair of same ori- 
ented spins, with probability 'Wo^ 

• Branching of a kink: 

can be generated by two overlapping spin-flips 
(with the probability wf ) such that it causes a spin- 
exchange. 

There are two independent parameters related to the 
original parametrization of Glauber: 



w, = r(i - 6)/2 

Wo = r(i + ~5) . 



(2) 
(3) 



In the present study we investigated the F = 0.35 and 
r = 2 cases. The corresponding PC critical points, 
without external field, are located at 5c = —0.535 and 
5c — —0.416 respectively. 

A single Monte Carlo step (MCS) consists of updat- 
ing all sites at once as described in [43] (throughout the 
paper the time is measured by MCS). The simulations 
were reformed on L = 40000 sized lattices, with periodic 
boundary conditions up to tmax = 10® MCS. The change 
of critical point 5c{h) has been determined for several 
values of external field. The transition probabilities are 
modified in the presence of an external magnetic field H 
as: 



..h _ 



w: = Wi{l - hsi), 



w;. = Wo{l - hsi), 
H . 



th{ 



kT' 



(4) 
(5) 

(6) 



Fig-fflshows the phase diagram in the (ft,, 5) plane, which 
is similar to the one we obtained in [24] by simulations 
and cluster mean-field technique. We have applied only 
random initial state simulations to find the points of the 
line of phase transition exhibitingkink density decay ex- 
ponent a — 0.1595 of the DP [3|. A clear power-law 
scaling of the critical point shift A5 as the function of h 
can be fitted with the form ^ in the region Q < h < 0.01 
with l/(j> = 0.52(3) crossover exponent. The error-bars 
on Fig. [1] come from our numerical estimates for the crit- 
ical point shifts and the least squares error estimate fit- 
ting procedure. This value is close to the early numerical 
estimates: = 2.1(1) by and = 2.24(10) by 
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FIG. 1: (Color online) Critical point shift in the NEKIMCA 
model as a result of the external h field for F = 2. The solid 
line is power-law fitting on the numerical data points. The 
insert shows the scaling of the shift of 5 for F = 2 (bullets) 
and for F — 0.35 (stars) on a log-log scale. 
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FIG. 2: (Color online) Critical point shift of the Grassberger's 
A SCA model to DP (bullets). The lines show power-law 
type fitting on the numerical data. The insert shows the local 
slopes of the crossover exponent defined as Q. 



III. CROSSOVER OF THE GRASSBERGER-A 
SCA MODEL 

The PC conserving SCA by Grassberger is realized by 
the following range- 1 update rule (we show the config- 
urations at t — 1 and the probability p of getting '1' at 
time t): 

t-1: 100 001 101 110 Oil 111 000 010 
t: 1 1 1-p 1-p 1 

The time evolution pattern in 1 -I- 1 dimension, for small 
p evolves towards a stripe-like ordered steady state (with 
double degeneration), while for p > pc — 0.1245(5) the 
kinks (the '00' and '11' pairs) survive. According to the 
classification of Wolfram ^] for p = we have the Rule- 
94, class 1 CA, while the p = 1 limit is the chaotic Rule-22 
deterministic CA. This model has been investigated from 
damage spreading point of view by us (5l| . 

Now we extend this model in such a way that a Z2 
symmetry breaking occurs in the absorbing state. This 
can be achieved by favoring one of the absorbing phases 
('10101010' or '01010101') shifted by a single site. There- 
fore we modified the transition rates 



t-1: 100 001 101 110 
t: 1 1 1-p' 



Oil 111 000 010 
1-p' 1 



such that p' — p + w a,t odd and p' = p ~ w at even sites. 
The simulations were run on L = 10^ sized rings up to 
tmax = 10^ MCS. The location of the critical point p{w) 
is shown on Fig. O To get more precise estimates for the 
crossover exponent we determined the local slopes of the 
data points: 



Inrcjwj) - Inre(wi-i) 
ln{wi) - ln(wi_i) 



(7) 



This is plotted in the insert of Figl2] and in the w ^ 
limit one can read-off the l/(j> — 0.53(2) linear extrapo- 
lation value (with least squares error estimates) in agree- 
ment with the value for NEKIMCA-|-h. As the figure 
shows the correction to scaling can not be neglected. 



IV. CROSSOVER IN A BARWE MODEL 

To study the PC to DP crossover in an other way we 
take a siinple version of a BARWe particle model intro- 
duced in [Bj] and generalized in (ZAMb model) to al- 
low phase transition at finite branching rate. This model 
is defined on the one dimensional lattice as follows. An 
occupied site is chosen randomly and is tried for diffu- 
sion, with probability D, or branching with probability 
a = 1 — D); while the time is increased by 1/N, where 
N is the number of occupied (active) sites. In a diffusion 
step the particle jumps to its randomly chosen nearest 
neighbor sites. If the site is occupied both particles are 
annihilated with probability r. On the other hand the 
jump is rejected with probability 1 — r. The branching 
process involves the creation of two new particles around 
the neighborhood. If either, or both neighboring sites 
is previously occupied the target site(s) become empty 
with probability r. Otherwise the lattice remains unal- 
tered with probability 1 — r. 

To induce a crossover to DP we introduced a sponta- 
neous particle removal ^ ^ or a coagulation AA A 
with a small probability w to the above reactions. The 
simulations were done on lattices of size L = 10^ with 
periodic boundary conditions for diffusion rate I? = 0.2 
(the convergence of GMF approximations was found to 
be rather good for this diffusion rate in (E^l)- For 
w = we used the critical point value determined in 
[5^ : Tc = 0.562(1). As we increase w the critical point 
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rc{D) shifts as shown on Figure [3l A simple power-law 
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FIG. 3: (Color online) Critical point shift in the ZAMb model 
to DP as the result of AA A (bullets) and A ^ (boxes) 
parity conservation breaking. Dashed lines are power-law 
type fitting to numerical data. The insert shows the corre- 
sponding effective exponents defined as ([7]). 

fitting for the data in the range w € (0,0.1) of the form 
[T] resulted in Arc = O.lSOui"'^*^. A more precise esti- 
mate, which takes the corrections to scaling into account 
can be obtained by calculating the local slopes ([7]). By 
plotting the effective exponent as the function of w a, lin- 
ear extrapolation resulted in l/cj) = 0.205(5) both for the 
AA A and the A ^ parity breaking reactions 

V. GMF+CAM CALCULATIONS 



where the site variables take the values: = 0, A. Dur- 
ing the solution of these equations one estimates larger 
than N sized block probabilities by the maximum overlap 
approximation: 

„ , ^ Pn{si, ...SN)PNis2, ■■■SN+l) 

. . (9) 

Taking into account the spatial symmetries and the con- 
servation of probability for the maximal, N = 9 approx- 
imation of this work we had to find the solution of a set 
of nonlinear equations of 272 independent variables. 

The steady state solutions of the TV-cluster approxi- 
mations for the NEKIMCA (1 < TV < 8) and the ZAMb 
model (1 < A^ < 9) have been determined and the corre- 
sponding densities are calculated numerically. The phase 
transition points are obtained for several values of the 
crossover parameter in both cases. In case of the ZAMb 
model (at D = 0.2) they are plotted on Fig[4] As one 
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The simulation results were complemented by ana- 
lytical cluster mean-field approximations and coherent 
anomaly extrapolations. The generalized (cluster) mean- 
field method (GMF) suggested for SCA by H and for 
dynamical RD models by [5^ has been shown to be 
successful approach for exploring the phase diagram of 
nonequilibrium models (see for example (sgI . [57| and 
the references in Q). It's extension with the coherent 
anomaly (CAM) extrapolation enables one to extract 
the true scaling behavior. 

In GMF we set up equations for the steady state of 
the system based on A^-point block probabilities. Corre- 
lations with a range greater than A^ are neglected. By 
increasing A^ from 1 (traditional mean-field) step by step 
we take into account more and more correlations and get 
better approximations. One can set up a master equa- 
tions for the Pn block probabilities as 



^PN{{S^}) 

dt 



(8) 



Note that the sign of correction to scaling is different for the two 



FIG. 4: (Color onUne) Critical point shift of the ZAMb model 
to DP (bullets) obtained by the GMF method for various 
cluster sizes (A = 1, 2, 9 top to bottom). 

can see for A'' = 1 (site mean-field) the critical point 
rc{N,w) exhibits linear relation without corrections to 
scaling. This also true for the ZAMB-|-sink case. This 
means the l/(t>MF = 1 mean- field value for both kinds of 
crossover. 

For A > 1 the rc{N,w) curves pick up w^^^ type of 
corrections to scaling (the leading order singularity is ex- 
pected to remain mean- field like). Assuming a scaling 
with the correction form 



r^N, w) = a{N)w + h{N)w^/'^ 



(10) 



crossover versions 



which is natural if we swap the axes Tc and w, one can de- 
termine the amplitudes {a{N)) of the leading order term. 
Following the envelope scaling hypothesis of the Co- 
herent Anomaly Method (CAM) extrapolation can be 
performed on the amplitude data in the N ^ oo limit. 

According to CAM the amplitudes a{N) of the cluster 
mean-field singularities in the leading order scale as 

\aiN)\oc\r,{N)-r,\'/'^-'/^-- (11) 
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FIG. 5: (Color online) CAM extrapolation results for the 
anomaly amplitudes a{N). Bullets: NEKIMCA+h (3 < iV < 
8 right to left), boxes: ZAMb + sink (5 < < 9 right 
to left)). The numerical errors are smaller than the symbol 
sizes. The dashed lines show the expected critical behavior 
from simulations. 



FIG. 6: (Color online) Schematic mapping of the two differ- 
ent kinds of crossovers between the kink and the spin model. 
The top graph shows that the removal of a single kink cor- 
responds to flipping of all spins to the right from position 
"X". The bottom graph demonstrates that an external mag- 
netic field in the spin model causes repulsion/attraction force 
is among domain walls (horizontal arrows) and an effectively 
two-component reaction-diffusion system of kinks (particles) 
emerges. 



VI. CONCLUSIONS AND DISCUSSION 



allowing to estimate the 1 /(f> exponent of the true singu- 
lar behavior (Eq.([T|)). The a{N) amplitudes were deter- 
mined by a next leading order fitting form (jlO[) on the 
the rc{N, w) data in the neighborhood of rc{N, 0) for the 
ZAMb model. In case of the NEKIMCA the same proce- 
dure has been applied for the 6c{N,h) crossover critical 
point GMF results. 

The highest a{N) amplitudes with a CAM fitting form 
(which takes into account possible scaling corrections 

m) 



a(7V) = cA^ + dA^+i , (12) 



where Ac{N) = \6c~SciN)\ for NEKIMCA and Ac{N) = 
\rc - rciN)\ for the ZAMb are plotted on Fig. [5l In this 
form the nonuniversal fitting parameters are: c, d and the 
anomaly exponent is : x = By the simula- 

tions (apart from the constants c and d) we have an ex- 
pected behavior for ^ both for ZAMb and NEKIMCA. 
This is plotted on Fig. [5] by the dashed lines. As one can 
see our GMF data can be fitted with those lines, (with 
the values: c = 2, d = 4.6 for ZAMb and c = 0.21, d = 
for NEKIMCA) but the amphtudes for iV < 10 just start 
to converge towards the asymptotic scaling curves. This 
suggests that larger A^-cluster approximations should be 
determined. However the numerical instability of root 
finding of the GMF method (within the space of more 
than 500 variable, nonlinear system) prevented us to go 
further. Still the different behavior for NEKIMCA and 
ZAMb can be justified by this figure. 



We have performed simulations and GMF-I-CAM ap- 
proximations for various versions of the PC class model 
crossovers towards the " DP class" . By determining the 
crossover exponents we have found and outstanding dif- 
ference between the effect of a symmetry breaking field 
(1/0 = 0.53(2)) and the parity conservation breaking 
(1/0 = 0.205(5)). Although the NEKIMCA and the 
BARWe are dual models, the two types of crossover op- 
erators can't be mapped onto each other. Furthermore 
they can't be mapped onto a local one in the dual sys- 
tem. The A — > process in the spin language would 
mean flipping of all spins from a given site (see FiglH]). 
Such transformation can't even be done in finite systems 
with periodic boundary conditions. 

A Z2 symmetry breaking in the BARWe model favors 
the (up) or '— ' (down) oriented spins depending its 
sign. In case of a '-)-' preference (see FiglHl) this means 
that '-|-'domains broaden on the expense of '— ' domains, 
hence odd-even kinks (the corresponding particle pairs 
in BARWe) attract and even-odd pairs repel each other 
(shown by the horizontal arrows). Therefore that parti- 
cle model becomes an effectively two-component, parity 
conserving one, and the sufficient conditions of the DP 
hypothesis do not hold. Still one can see a DP type of 
decay in the global order parameter because calling a 
domain a macroscopic particle 'X', an effective DP pro- 
cess: X ^ ^, X ^ 2X, 2X X describes it's dynamics. 

However there are certain operators which do not ex- 
hibit DP type behavior. For example in case of cluster 
spreading it is easy to see that there are two sectors in 
this model. An odd parity one, with f3' = Q final survival 
probability defined as 

Poo(x\p-pY (13) 

and an even parity one, with normal, j3' = (3 DP expo- 
nent. As a consequence the hyper-scaling relation con- 
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necting the cluster spreading exponents and the rapidity 
reversal symmetry of DP is not satisfied for this critical 
behavior. The universality class of the DP is split into 
sub-classes. The subclass of models with BARWe dy- 
namics with broken symmetry (let's call it DP-2) is 
different from that of the ordinary 1+1 dimensional DP. 
Since the difference is manifested in the cluster spreading 
behavior, which is a consequence of a special initial con- 
dition, similarly to the terminology used in equilibrium 
models with different surface classes we do not claim the 
splitting of the DP class itself. 

This sub-class behavior is similar to the one, which is 
observed in the pair contact process, where the frozen 
monomers cause different cluster behavior from that of 
the simple contact process. Also that kind of DP-2 be- 
havior should arise, when one introduces spin-anisotropy 
in the NEKIM [H, [G^I and the effective two-component 
BARWe model exhibits DP type of phase transition for 
finite branching rate in the global order parameter, but 
different scaling behavior occurs in terms the cluster 
spreading. Furthermore this model exhibits a re-entrant 
phase diagram and at zero branching rate one finds an- 
other phase transition belonging to the two-component 
BARWe model [3l|, but the cluster spreading behavior, 
which is sensitive to the spin anisotropy is different again. 
Therefore such sub-classes are not at all rare among low- 
dimensional reaction-diffusion models. 

We would like to point out that by using an alterna- 
tive definition for universality classes a different interpre- 
tation can also be given for the same numerical results 
presented here. Some authors define a university class 
by the field-theoretic action of the model (without spec- 
ifying the fixed point) instead of the models exhibiting 
the same set of known critical exponents as we do. Since 
the action of BARWe particles |3l| and the action of the 
NEKIM spin model [34I (called GV model) do not agree, 
they should belong to different universality classes, which 
"intersect" in one dimension only [2]. According to this 
picture the applied symmetry breaking to the PC and GV 
class models (.Z2 breaking in GV and particle removal in 
BARWe) should end up in DP class somehow. 

However we think that our definition of universality 
class, and therefore our interpretation, is more precise. It 

^ Note that this is not a finite size effect 



can describe the critical behavior of (particle) system for 
which no proper field theoretical action has been found. 
Even within a field theory multiple fixed point solutions 
(corresponding to multiple critical points) can exist and 
the relevancy of terms affecting the stability of a fixed 
point is not clear in many cases. 

Furthermore a coarse grained field theory may not cap- 
ture all scaling details of a particle model. In a discrete 
particle model diffusive annihilating particles can die out 
within finite time (even in an infinite sized system) due 
to the recurrence relation of random walks in one dimen- 
sion ^ . The particle survival probability decays asymp- 
totically as Ps(0 ^ ^-iid the final survival probability 
scales as (jl3|) . As a consequence a thermodynamic limit, 
in which the particle density are kept finite can't be es- 
tablished. 

On the other hand in field theories of continuous vari- 
ables the survival probability is always unity and it had 
been unclear if the exponent 5 and /?' could be defined 
sensibly. To overcome this discrepancy a reinterpretation 
of survival probability was proposed in [6l[ . This alterna- 
tive definition of (5 in field theory resulted in correct scal- 
ing exponents and hyper-scaling relations corresponding 
to symmetries. 

Since in our case the "thermodynamic limit" can't be 
restricted to infinite particle number the cluster spread- 
ing behavior starting from finite number of particles is 
relevant from the universality point of view as in many 
papers investigating absorbi ng p hase transition via this 
approach (see [E H, i, i, [uf [H, [13. The cluster 
spreading behavior (the finite survival probability, pair- 
connectedness functions, avalanche distributions ... etc. 
) are sensitive to the initial parity of particles, hence 
strictly speaking the DP-2 subclass is not identical to 
the DP class in 1+1 dimension. 
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